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CKM Physics from Lattice QCD 

Paul Mackenzie 

Fermi lab, P.O. Box 500, Batavia, IL, USA 



I discuss the lattice calculations relevant to recent advances in CKM phenomenology, focusing on those relevant 
to new experimental results reported at this workshop. 
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1. Introduction 

Lattice QCD calculations make possible some of 
the most important results in CKM phenomenology. 
For example, with lattice QCD, new constraints on 
the CKM matrix can be obtained from the two most 
important new experimental results reported at this 
conference: the observation of the purely leptonic de- 
cay of the B meson by Belle 1, 2j, and the measure- 
ment of Bs'B's mixing by CDF 3] and DO 4]. These 
quantities are examples of golden quantities for lattice 
QCD, processes with stable particles and a maximum 
of one hadron in the lattice box at a time. Other 
such quantities relevant to CKM phenomenology in- 
clude other meson leptonic and semileptonic decays 
and MM mixings. Good prototype calculations exist 
in the quenched approximation (ignoring light quark 
loops), and many unquenched. We know of no im- 
pediment to making them more and more accurate as 
methods and computers continue to improve. 

There are several different methods for putting 
quarks on a lattice that have widely differing virtues 
and defects. Wilson and clover (or improved Wilson) 
fermions have been the most widely used. They em- 
ploy a strong breaking of chiral symmetry at the lat- 
tice spacing scale for technical reasons (to solve the 
fermion doubling problem). This has lead to diffi- 
culties in reducing the light quark masses to close to 
their physical values. Staggered fermions have a chi- 
ral symmetry. This makes possible closer approach 
to the physical, light quark mass limit, smaller uncer- 
tainties from chiral extrapolation and statistics than 
is possible with other methods. However, they also 
have residual fermion doubling, which in unquenched 
calculations is cured by taking a root of the fermion 
determinant. This procedure has been shown to make 
sense in every calculational test to it has been put, 
but more work is on-going to investigate it. Two new 
related methods, domain wall fermions and overlap 
fermions, suffer from neither of these defects, but are 
computationally much more expensive. Practical cal- 
culations with these last two methods are in their in- 
fancy. There are therefore multiple opinions among 
lattice theorists about the optimal fermion method 
for phenomenological calculations at present. Several 
methods may have to be pushed to much higher ac- 
curacy to test that all methods give the same results 
and resolve the question. 
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Figure 1: The light quark mass dependence of 
fss \J I {/b V Mb) with improved staggered fermions 
(HPQCD, nf — 2 + 1, black symbols) and clover fermions 
(JOQCD, n/ = 2, red symbols). Black and red lines are 
extrapolations using the two sets of data. The vertical 
dashed line is the physical light quark mass. ^ 



In this talk, I will discuss mainly lattice results of 
direct interest to this workshop. For recent more com- 
plete reviews, see jS»i2ei2iJ- 



2. Meson Decay Constants 

The practical trade-offs between Wilson and stag- 
gered fermions are illustrated in Fig. ^ It shows 
/b^ \/MBj{fB^MB) as a function of light quark 
mass. The vertical dashed line shows the physical 
quark mass. The Wilson data 8] are shown in red, 
and the staggered data 9] in blue. The two methods 
agree where their data overlap, at mi > but 
curvature at light quark masses yields a significantly 
higher value at the physical quark mass than is ob- 
tained from a naive extrapolation of the Wilson data 
from the region mi > ms/2. The staggered data yield 
the result 



/b = 0.216(9)(19)(4)(6) GeV. 



(1) 



The uncertainties in and /b^, other than those 
arising from chiral extrapolation and statistics largely, 
cancel out in the ratio. HPQCD finds, with staggered 
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fermions ^ in blue., 

fBjfB = 1.20(3)(1), 



(2) 



where the first error is from chiral extrapolation and 
statistics and the second is from everything else. This 
can be compared with the result that JLQCD found 
earlier, chirally extrapolating from a much larger 
quark mass: JeJIb = 1.13(3)( ^2')- 



For /b,, HPQCD quotes 



/b^ = 0.260(7)(26)(9) GeV. 



(3) 



The uncertainties in /bJ/b and /b^ are relatively 
independent of each other, since the uncertainties in 
the former are dominated by statistics and chiral ex- 
trapolation, while the uncertainties in the latter are 
dominated by everything else. /bJIb is more inde- 
pendent of /b^ than of /b- 

D and Dg meson decay constants are also of great 
interest because of the improvements to the exper- 
imental numbers coming from CLEO-c. The issues 
for these calculations are very similar to those for the 
B mesons. With staggered fermions and 2 -h 1 un- 
quenched light quark flavors (that is, degenerate up 
and down quarks, and a nondegenerate strange), Fer- 
milab/MILC obtain 0| 

fn = 201(03)stat(17)sys MeV (4) 
fD^ = 249(03)stat(16)sys McV. (5) 

With clover fermions and 2 light flavors, CP-PACS 
obtains 

+20 

fo = 202(12),tat(_23)sysMeV (6) 



+07 

fo^ = 238(ll)stat( _27)sysM^V- (7) 



This may be compared with the recent result from 
CLEO-c 111 



fo = 223(17)(3) MeV, 



(8) 



assuming the canonical value of Vud- Alternatively, 
the comparison could be used to obtain an indepen- 
dent determination of Vud- 



3. Semileptonic Decays 

The shape of the form factor for the decay D ^ Klv 
was obtained in 2004 by the Fermilab/MILC collab- 
oration using unquenched staggered fermions It 
has been confirmed by several experiments with in- 
creasing accuracy, most recently by Belle Ts^ and 
BaBar. (See Fig. [21) CLEO-c is expected to deter- 
mine this shape with even greater precision, yielding 
an even more stringent test. 




4 experiment [Belle, hep-ex/05 1 0003 ] 
lattice QCD [Fermilab/MILC, hep-phy0408306] 
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Figure 2: The predicted shape of the form factor for 
D ^ Klv has been tested with increasing accuracy by 
experiment, most recently by Belle [l^. 



The decay B ^ -kIv is measure of Vuh^ and pro- 
vides a measurement of the height of the unitarity 
triangle that is competitive with h ^ u inclusive de- 
cays. For B ttIu^ the shape of the form factors is 
a critical issue, even for extracting the CKM matrix 
element. The reason is that the uncertainties in the 
form factors are highly q'^ dependent for both theory 
and experiment. In lattice calculations, predictions 
do not yet exist in the high recoil region. 

It has long been knows that analyticity and unitar- 
ity can be used to constrain the shape of form factors. 
Arnesen et al. have recently applied these ideas 
to B decay. They show that when the variable q'^ 
is mapped into a well chosen new variable, unitar- 
ity constrains the coefficients in the z expansion of the 
form factors so that only five or six terms suffice to de- 
scribe the form factors to 1% accuracy. They propose 
to to use various theoretical methods to determine the 
coefficients, including lattice QCD. 

Becher and Hill have recently pointed out that 
heavy quark theory leads to the expectation that 
in the heavy quark limit, the coefficients fall as 
(Aqcd/Mh)"*^'^, an even tighter constraints on the 
form factors, jl^ They argue that this leads to ex- 
pectation that only two or three terms are required to 
describe the data to high accuracy, and show that the 
BaBar data so far conform to this expectaion. (See 
Fig. El) 

This expectation is also confirmed in lattice data re- 
analyzed the same way 15j. Therefore, in this expan- 
sion, to the accuracy of the present data, comparing 
theory and experiment means testing theory against 
experiment with the slope, comparing the overall nor- 
malizations to obtain Vub^ and searching for evidence 
of the curvature term. 
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and 




t (GeV ) 




(10) 



For his best estimate of the combined matrix ele- 
ments, he combined these with the HPQCD, staggered 
fermion results for the decay constants already dis- 
cussed, and obtained 



fBVBB = 244(26) MeV 



(11) 



and 



^ = /B.//B\/4/^=1.210(+gJ). (12) 
Following this same procedure, one can obtain for the 



Ib^^Bb, = 294(33) MeV. 



(13) 



This quantity and ^ are the combination of BB mixing 
quantities with uncertainties that are most indepen- 
dent of each other. 

Lifetime differences require a different operator 
than the one for the mass difference. Unquenched 
calculations of matrix elements for this operator do 
not yet exist. The required operator has been given 
as part of the complete set of four-quark operators, in- 
cluding those arising in supersymmetric theories, and 
initial results have been reported in the quenched ap- 
proximation [r^ . 



Figure 3: Plotting the form facto shape for B ^ nlv decay 
versus a well chosen new variable z (instead of t = pro- 
duces a function that is linear to the accuracy of existing 
data, and can thus be describe with only two parameters. 



4. Meson Mixings 

Matrix elements for BB and BgBg mixing are con- 
ventionally parameterized in terms of the meson de- 
cay constants and bag parameters, as in the kaon sys- 
tem. What is actually calculated on the lattice, how- 
ever, are the combinations /bV^b and /bs^/B^. 
Since, these combinations are what is needed for phe- 
nomenology, these are what lattice theorists should be 
reporting. 

Unquenched results were reviewed at Lattice 2005 
by Okamoto [51] • Since unquenched results using stag- 
gered fermions had not yet appeared, he took results 
for bag parameters from JLQCD using two flavors of 
clover quarks. 



5. Effect on the p - plane 

5.1. Effects of Lattice Fermion Methods 

In this talk, I have emphasized staggered fermion 
results because they are for the most part fully un- 
quenched, with the right number of light quark fla- 
vors, and with light quarks masses much closer to their 
physical values. Other lattice theorists have a prefer- 
ence for other fermion methods, and the lattice com- 
munity is not yet settled on the optimum approach. 
Another widely used set of numbers for BB mixing pa- 
rameters comes from the CERN CKM study of 2003 
^IR] (based on the 2002 review of Lellouch ^]) 




235(33)(24) 
276(38) 



C = 1.18(4)(^^^). 



(14) 

(15) 
(16) 



B(mb)= 0.836(27)(_g2) 



(9) 



These were based on a set of mostly clover lat- 
tice results, most in the quenched approximation but 
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Figure 4: Constraints on the p — rj plane derived from B meson leptonic decay and BB mixing. The bounds in the left 
figure employ clover fermions, those on the right employ staggered. 



with some two flavor results used to extrapolate to 
the three-flavor theory. The results are compatible 
with the staggered results already discussed, but with 
larger uncertainties. 

The effects of the smaller uncertainties of staggered 
fermions were shown by J. Charles at this workshop 
[19]. Fig. SI shows the bounds on the p — rf plane de- 
rived solely from B meson leptonic decay and BB mix- 
ing. The bounds in the left figure employ the clover 
fermion results just discussed, those on the right em- 
ploy the staggered fermion results discussed earlier. 



5.2. CKMf itter and UTFit 



There are two groups producing widely cited global 
fits to CKM data: CKMfitter ^ and UTfit |2ll. 
The groups employ different statistical methods, and 
obtain somewhat different results. For example, for 
Amg predicted without incorporating the experimen- 
tal BgBg mixing results, UTfit obtains 21.5(2.6) ps~^, 
while CKMfitter reports 21.7(4 2) P^"^- ^his is puz- 
zling, since one would expect this quantity to be 
sensitive mainly to lattice uncertainties, and both 
groups take the lattice results of Ref. 0| as the 
starting point. One possible difference is that UT- 
fit lists fs^ \J Bb^ and ^ as inputs, while CKMfitter 

lists fB\/ Bb and ^. The latter combination contains 
highly correlated uncertainties that must be treated 
with care. However, a more significant difference 
seems to arise from statistical methods and treatments 
of combinations of uncertainties, rather than differ- 
ences of lattice inputs. These differences remain yet 
to be resolved. 



6. Outlook 

Lattice calculations are playing an essential role in 
enabling some of the most important results in the 
CKM experimental program. Most of the key lattice 
calculations for CKM physics involve single-hadron 
processes with hadronically stable meson, which are 
among the most solid current lattice calculations. Sig- 
nificant issues are currently being worked through, for 
example, the best fermion method for phenomeno- 
logical calculations, the best way to incorporate lat- 
tice results into Standard Model global fits, and the 
best way compare theory and experiment in semilep- 
tonic decays. Lattice methods are currently in a state 
of productive ferment, with several different methods 
for unquenched lattice fermions under active investi- 
gation by various groups. Lots of progress is being 
made in algorithms for these various methods g22j . 
The computing power being applied to lattice phe- 
nomenology is rising exponentially 23]. There are ex- 
cellent prospects for further progress in lattice CKM 
phenomenology. 
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